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Abstract: We use the techniques of effective field theory in an expanding universe 
to examine the effect of choosing an excited inflationary initial state built over the 
Bunch-Davies state on the CMB bi-spectrum. We find that even for Hadamard states, 
there are unexpected enhancements in the bi-spectrum for certain configurations in 
momentum space due to interactions of modes in the early stages of inflation. These 
enhancements can be parametrically larger than the standard ones and are potentially 
observable in future data. These initial state effects have a characteristic signature in 
/-space which distinguishes them from the usual contributions, with the enhancement 
being most pronounced for configurations corresponding to flattened triangles for which 
two momenta are collinear. 
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1. Introduction 

Inflation has become the dominant paradigm for the study of the early universe. Cur- 
rent data from the CMB power spectrum [1], as well as large scale structure [2] are 
certainly consistent with the assumption of an inflationary phase in the early universe, 
while the anti-correlation between the TT and TE power spectra on superhorizon scales 
should soon become strong enough statistically to serve as a "smoking gun" for infla- 
tion. Furthermore, the case for inflation will certainly be strengthened by the discovery 
of the B-mode polarization of the CMB. 

These cosmological observations provide a powerful lever arm that allows us access 
to physics at energy scales far beyond those that will be probed by accelerators in the 
foreseeable future. Already, we can use measurements of the spectral index n s and 
bounds on the ratio of tensor to scalar fluctuations r to place restrictions on the form 
of the inflaton scalar potential, at least within the field range corresponding to the 
observable 10 or so e- folds of inflation [1]. 

The inflaton potential is only part of the story, however. The inflaton should al- 
most certainly be viewed as an effective degree of freedom, perhaps arising from a 
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higher dimensional theory as the separation mode between a brane and an antibrane 
[3], or as one of the many axions that appears in string theory [4], amongst many other 
possibilities. As such, the inflaton effective action will contain not only renormaliz- 
able terms, but also so-called irrelevant operators which are suppressed by powers of 
the cutoff scale M. At energies higher than M this effective theory breaks down and 
the inflaton should not be viewed as the appropriate degree of freedom. The effect 
of these operators needs to be taken into account when computing the various infla- 
tionary observables, and several works have considered these corrections [5]. Whilst 
these corrections might a priori have been expected to be small, they can give rise to 
corrections to the CMB bi-spectrum which are larger than the standard result as shown 
by Creminelli [6]. This arises because the standard contribution to the bi-spectrum is 
already highly suppressed by the slow roll parameters [7, 8]. 

Equally as important as including higher order terms in the effective action is a 
proper treatment of the quantum state of the fluctuations of the inflaton about its zero 
mode. The standard calculations of the power spectrum are predicated on a particular 
choice of quantum state, the so-called Bunch-Davies (BD) state [9]. There has been 
a great deal of work in recent years focussed on understanding what might constitute 
reasonable modifications to this quantum state and how the data we already have on 
inflationary observables from the power spectrum might constrain these modifications 
[10]. 

In this work we ask the question: can we probe the nature of the initial quantum 
state of the inflaton through the non-gaussianities produced via inflaton fluctuations? 
We will argue that statistics that probe this non-gaussianity are amazingly sensitive to 
the nature of both the interactions of the inflaton and more importantly, to its initial 
state. In fact, in some ways, higher correlation functions are more sensitive probes 
of initial state fluctuations than the power spectrum. The reason for this is that the 
power spectrum is only sensitive to the interactions of the inflaton through loops which 
are highly suppressed [If, f3] (see also [12] for slightly different conclusions). On the 
other hand, non-gaussianities can probe the inflaton's interactions directly at tree-level, 
which can be significant if the inflaton does not start out in the BD vacuum state. 

The simplest correlation function that probes the non-gaussian nature of inflaton 
statistics is the three point correlation function, otherwise known as the bi-spectrum, of 
the fluctuations. The seminal calculation of the bi-spectrum was done by Maldacena [7] 
(see also [8]) who expanded the action for a minimally coupled inflaton scalar coupled 
to gravity to third order in the gauge invariant curvature fluctuation variable (. This 
is the relevant quantity since it remains constant outside the horizon, re-emerging as a 
non-gaussian contribution to perturbations on scales relevant to present day cosmology. 
Subsequent calculations have refined estimates of the bi-spectrum in the standard case, 
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extended them to other inflationary models [14], and have also considered the effect of 
irrelevant operators [6]. 

Let us briefly sketch the details of the calculations. In a cosmological setting, we 
need to consider the time evolution of equal time correlation functions. Thus we choose 
an initial state at a (conformal) time r) , which we take be at or near the onset of the 
inflationary phase. There is considerable freedom in how to choose this state and we 
will discuss this issue further below. To compute the time evolution, we work in the 
interaction picture where the state evolves according to 

;*^>>=W(,)>, (i.i) 

where Hi is the interaction Hamiltonian. This has the formal solution given an initial 
state \ip{j]o)) at a (finite) time r\ = i] : 

mr ] )) = Te- t Vo H < WW mr ]o )}. (1.2) 

We require that our effective theory be valid at i] = i] . In particular, as will be discussed 
in Sec. 2, this means that states are not excited at momentum scales greater than the 
cutoff. The initial state can be described by giving all of the correlation functions of 
our dynamical degrees of freedom, e.g. the comoving curvature perturbation C( 7 ?)- 111 
particular the tree level contribution to the equal time three point function is given by 

^(v)\c^kM(k 3 (v)mv)) = mvomsvx^kM^o)} (i.3) 

- i /%V^(^0)|[C£ 1 (^)C£ 2 (^)C4(^)>^(V)]l^(^0)) + O(^ 2 ). 

Thus the total three point correlator is a sum of a contribution from any initial non- 
gaussianity present in the state at the beginning of inflation, evolved forward in time 
with the free Hamiltonian, and the contribution that occurs from interactions that 
take place between the beginning of inflation and the time of observation. The former 
contribution has to be treated as arbitrary until we have a better understanding of 
the pre-inflationary stage. At best, we can place bounds on its contribution based on 
ensuring that its backreaction is small. The second contribution, on the other hand, 
depends explicitly on the interactions of the theory as well as the initial state. This 
contribution will be the focus of our investigations in this work. 

In practice we shall assume that the initial state is gaussian, so that all the corre- 
lation functions of operators in this state will be completely specified by the two point 
Wightman function (non time-ordered in- in expectation value), which in momentum 
space takes the form 

mVo)\(j;Sv)(j;M)\^o)) = (27rf5%k 1 + k 2 )G> i ( V , V ') = ^'(h + k^MV^W), 

(1.4) 
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Figure 1: Momentum triangles which dominate the bi-spectrum: 1. Equilateral triangles 
arising from Hubble- horizon crossing. 2. Squeezed or local triangles arising from super-Hubble 
evolution. 3. Flattened triangles arising from intial state effects. 

where Vh are properly normalized solutions of the linear equations of motion for ( 
in momentum space. In Sec. 2 we discuss the conditions required of this two-point 
function so that we can reliably trust our calculation. 

In the standard calculation using the BD vacuum, the dominant contribution to the 
three point function, which measures directly three particle interactions, comes from 
when the modes cross the horizon. The intuitive reason for this is that at sub-horizon 
scales the BD vacuum corresponds to a state of no-particles 1 (where the particles are 
inflaton quanta). Thus there are no particles to interact and so no contribution to the 
three point function. As the modes cross the horizon, the WKB approximation breaks 
down which is tantamount to the statement that particles are created. These particles 
can then undergo interactions which contribute to the three point function. Once the 
modes are well outside the horizon, (, properly defined, gets frozen in nonlinearly i.e. it 
is conserved and further interactions become irrelevant. In curvaton models [15] there 
can be an additional contribution from super-horizon scales because ( is not necessarily 
conserved [16]. In figure 1 we show the momentum triangles for which the three point 
function is maximized for these two types of effects (triangles 1 and 2). 

The main point of this paper is that this situation changes dramatically if the initial 
state is not BD. In this case there are particles present initially which can undergo 

1 It is true that for a static observer the BD vacuum appears as a thermal bath of particles. However 
here we are using the adiabatic or WKB definition of particles appropriate to the flat slicing of de 
Sitter. This is the definition that is actually most useful in the context of calculations. 
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interactions. Furthermore the interactions of the inflaton are necessarily stronger at 
the beginning of inflation than at the end, and so we get a second contribution to the 
non-gaussianity coming from interactions in the early stages of inflation. These effects 
dominate the bi-spectrum for the flattened triangles (defined as those for which two 
momentum are collinear) of type 3 in figure 1 (these were called 'folded' in Ref. [17] 
where a similar study was performed). To see why interactions get stronger in the 
past, let us consider the prototypical example of an interacting scalar field <fi on an 
FRW geometry ds 2 = a 2 (r])(—dr] 2 + dx 2 ) with Lagrangian expressed in conformal time 

■>> 2 - 1 - 1 (>)')" + ■■■)■ 

n=2 m=2 V 7 / 

(1.5) 

where M is the cutoff scale and the ellipsis corresponds to additional higher derivative 
interactions expected on usual grounds. Upon rescaling to the canonically normalized 
'comoving' field x — a( P we obtain 2 

s = fd v [ d 3 x \--{d x f + —x 2 -Y , A " A x n 

J 'J I 2 V AJ 2a ^ n!a"- 4 M™- 4A 

- £ M t- lal L ((Sx) 2 + - l a -xx) m + ...). (1.6) 

Note that the scaling with a is the same as that with the cutoff M. At subhorizon scales 
the terms suppressed by a' /a and a" /a are irrelevant. It is apparent that higher order 
potential interactions with n > 4 and higher order kinetic interactions with m > 2 
necessarily grow in the past as a — > 0. In other words the irrelevant operators become 
more relevant as modes become blue-shifted. In this case these terms will not contribute 
to the three point function at tree level (the first interesting contribution is to the 4-pt 
from the kinetic interactions) but as we shall see later when expanding around the 
background inflaton solution 4>o{rj), we easily generate contributions to the three point 
function that scale as positive powers oil/ a. Somewhat to our surprise then, it appears 
that the theory with higher derivative terms gives rise to a greater enhancement than 
the theory without them; the standard expectation is that irrelevant operators would 
be subdominant for low energy observables such as the CMB bi-spectrum 3 . 

The upshot of our calculation is that the bi-spectrum is at least as good a probe 
of the initial state of inflaton fluctuations as the power spectrum. The effects in both 

2 In the full gauge invariant calculation this is the familiar Mukhanov variable v. 

3 This is a reflection of the fact that higher derivative operators are irrelevant in the IR. However, 
the scales that are in the IR today, relative to the cutoff scale M, were in the UV, or at least less in 
the IR at earlier times due to the expansion of the Universe. 
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cases will depend on the Bogoliubov (3k coefficient (see below) but there is no extra 
enhancement for the case of the power spectrum. We also see that a great deal of 
physics can be missed by not including the effect of irrelevant operators. Earlier works 
[18] have also considered the effects of initial state effects but not directly there impli- 
cations for interactions. More recently in Ref. [19], similar but less explicit conclusions 
were made about the non-gaussianities. However, our approach to considering the typ- 
ical corrections and in particular the magnitude of the backreaction [20] is somewhat 
different. The general features of the effects considered here including the precise shape 
dependence were pointed out in Ref. [17], although explicit estimates of the possible 
magnitude of these effects were not given. 

In Sec. 2 we discuss how effective field theories are to be construed in the context of 
an expanding universe, and give a preliminary discussion of the importance of bounding 
backreaction. We then turn to our calculations of the bi-spectrum for the remormaliz- 
able and higher derivative interactions of the form described above in Sec. 3. We will 
see that for certain triangles in momentum space, there are additional enhancements 
as a direct consequence of the inflaton interactions taking place at the beginning of 
inflation. In Sec. 4 we then show that these effects have sufficient measure, so that the 
enhancement survives even after the bi-spectrum is converted to spherical harmonic 
space. This is where the higher derivative operators triumph over the renormalizable 
ones, even when the effects in both theories are calculated with the same excited states. 
In Sec. 5 we discuss the form of the enhancements from higher irrelevant operators and 
finally we conclude in Sec. 6. 

2. Effective Field Theory in an Expanding Universe 

The usual approach to effective field theories [21] consists of the following steps. First 
identify applicable regime of energies/temperatures in which we want to do physics. 
Next, decide what the relevant degrees of freedom in this regime should be, as well 
as what symmetries their dynamics should obey. Finally write down the most general 
lagrangian in terms of these degrees of freedom that incorporates the required symme- 
tries. Once this is done, observables of the theory can be calculated in this regime. In 
this calculation, a determination of how accurate the results must be (to match exper- 
imental data, say) tells us which terms and in particular, which irrelevant operators to 
keep in our action. 

While this approach is sufficient for standard particle physics observables, it is de- 
ficient when it comes to calculating in an expanding universe. In this situation, we are 
not calculating S-matrix elements (which may not even exist in some cosmologies, such 
as de Sitter [22]) which depend on boundary values at asymptotic times. Instead, we 
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need to compute the time evolution of correlation functions which requires us to com- 
pute so-called "in-in" matrix elements [23], and which incorporates the time evolution 
of the states, as well as the operators. In essence, we are solving a initial- value problem 
rather than a boundary one. 

The notion of an effective theory presupposes an energy scale M such that we 
are constrained to do physics only below this scale. The new twist arising from the 
universes' expansion is the fact that it induces a redshift of energies so that scales 
that were once larger than M will eventually become part of the low energy (< M) 
spectrum. In fact, we can turn this around to make the following statement. Consider 
a physical momentum scale k corresponding to a length scale on the CMB sky today. 
Then, at a (conformal) time r/ such that k/a(rj ) ~ M we reach the limit of validity 
of the effective theory for this scale. If we want our effective theory to be valid for all 
of the relevant scales in the CMB sky, we must then impose a limit on how far back 
in time we can trust this theory; we will let r] denote the earliest time at which the 
effective field theory can be trusted, and we shall take this time to be the beginning 
of inflation, although more generally we only require it to be the time at which the 
physical scales of the CMB today are of the order of the cutoff. 

The scale M can also infiltrate the description of the initial state of the fluctuations. 
The standard way this state is chosen is by arguing that at short enough distances, the 
fluctuations behave as if they are in flat space, and their vacuum state would reflect 
this. In practice, the 77 — > —00 limit of the solutions to the mode equation is taken, and 
the linear combination of the solutions that approaches a positive energy plane wave in 
this limit is then chosen. The state picked by this procedure in the case of a de Sitter 
universe is the Bunch-Davies (BD) state [9] ; we will denote these modes via U k in the 
following. 

Given that at times earlier than t]q it may not be permissible to treat the inflaton as 
a well defined degree of freedom and even if we could, we would certainly not be privy 
to the relevant dynamics, it does not seem reasonable to use a state whose definition 
required going to arbitrarily short distances. In fact, using this state requires the radical 
assumption that the description of the physics in terms of the inflaton fluctuations as 
a free scalar field in an FRW background is a valid one at all scales. 

A more reasonable description of the state should depend more explicitly on the 
domain of validity of the theory. Thus we set initial conditions at i] and write 

V k (v) = a k U k (ri) + f3 k U* k ( V ), \a k \ 2 - \(3 k \ 2 = 1, (2.1) 

i.e. the mode V k is a Bogoliubov transform of the BD state. As is well known, this can 
be viewed as an excited state built upon the BD state; the number density of particles 
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of momentum k is \j3k\ • The Bogoliubov coefficients encode information about the 
initial conditions satisfied by the relative to the BD modes: if we specify that 

Vk(Vo) = -iu7kVk(Vo), Uk(Vo) = -iwfcZ4(%), (2.2) 

then we see that 

Pk _ u k - Wk ^ 2 3) 

Ctfc Uk + yJk 

How should we fix f3k (or equivalently, Wk)l In Ref. [24] Wk is described in terms 
of a Laurent expansion in ou^/M, and a new renormalization procedure has to be 
implemented to account for the spacelike boundary on which the initial conditions are 
specified. In this work we follow the more standard approach and demand that the 
state constructed from the modes Vk be Hadamard [25]. This fixes the short-distance 
behavior to be the usual one, namely that (3k must fall of faster that 1/k 2 , but it may 
otherwise be arbitrary. We implement this by noticing that since our description in 
terms of an effective theory breaks down at the scale M, we should not excite any 
modes with energies higher than this. Thus we demand that (3k — > for k > Ma(r)o). 
Note that this assumption shields us from the transplankian problem [26] since as the 
universe expands and transplankian modes redshift to cisplankian scales, they enter in 
their vacuum state (3k — 0, and so only the subsequent cisplanckian dynamics will be 
important. 

Backreaction 

We can put some bounds on how large the non-vanishing (3k s can be by considering the 
issue of backreaction. This was treated within the context of the standard renormal- 
ization procedure in Ref. [27] and in Ref. [20, 24, 28] within the context of boundary 
renormalization. The basic issue is whether the energy density coming from the "par- 
ticles" of the BD vacuum contained in the initial state might overwhelm that of the 
inflaton zero mode, and thus prevent inflation from occurring. 

This energy density can be estimated by taking the crude model (3 k ~ (3 oe - k2 /( Ma ( r )o)y 
and demanding that the energy density of the nearly massless quanta of the inflaton 
be less than M^H 2 : 

1 f d 3 k 2 «M 4 ]Q |2 ,,4 < M 2 tt2 . |fl ,2 < M pl (t) ^ 

P~- 4 J 7-^ Ifcl fc ~^)TlA>l M <M pl H =Hft| <^T- (2-4) 

More generally we can imagine any smooth fall-off as considered in Ref. [27], where a 
fall-off of the form \(3 k \ 2 ~ 0(k 4+s ) is assumed, and 5 > to ensure the Hadamard 
condition. This gives rise to similar results (up to factors of 0(1)), and so for the 
purposes of this article we shall consider the above crude model. 
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The stress-energy tensor also contains contributions from interactions, and these 
contribution are non-zero as soon as we generate some non-gaussianities. We shall 
perform more explicit calculations of these contributions in Sec. 3.3, but it is straight- 
forward to estimate the bounds on these effects. On dimensional grounds, as long as we 
remain in the regime of effective field theory, the largest the contribution to the energy 
density (expanded in powers of (3$) can be is of order |/3o| n M 4 for n > 1. For \(3 \ < 1 
we naively expect the terms linear in (3q to dominate. However, any term which is odd 
in Po will also contain some powers of e 2lkr] inside the momentum integrals. The same 
situation already arises in free theory where we get a contribution to the energy density 
of the form 

Ap = -1 /|^ fe *e- + ,,. (2.5) 

The crucial point is that at early times, i.e. large 77, the rapid oscillations of the expo- 
nential damp its contribution. For instance for the model (3k ~ p Q e~ k2 ^ Ma ^ Vo ^ 2 , this con- 
tribution behaves as g-V^a^o) 2 _ e -4M 2 a(r )0 ) 2 /(^«(»y)) 2 , At early times the exponential 
is negligible, and once the spacetime has inflated to the point that a(f])/a(i] Q ) ps M/H 
and the exponential is of 0(1), the 1/a 4 scaling suppresses the energy density to be 
of order H 4 and consequently negligible. Thus even with interactions included, the 
backreaction to the energy density can be no larger than \(3 \ 2 M 4 . This is borne out 
by the more detailed calculation in Sec. 3.3. 

We must also make sure that the slow roll conditions are not violated. Since 
H = -tH 2 = -^^(p + p) and H = 2er]'H 3 = ~^r~ (p - SH(p + p)), where 

mpl mpl 

e and rj' are implicitly defined slow roll parameters, then assuming Ap ~ \(3 \ 2 M 4 , 
Ap ~ |/3 | 2 M 4 and Ap ~ |/3 | 2 #M 4 we have the bounds 

IAI < V7^1, (2.6) 
IAI < V^'^f- (2-7) 

This conclusion, also reached in [28], differs strongly from the conclusions of [19, 20]. 
That author gives more pessimistic estimates based on the order (5k contributions which 
neglect the oscillating nature of the exponentials. It is also assumed that Ap ~ M Ap 
which again is not the case because of the oscillatory damping; rather we obtain Ap ~ 
HAp. 

We see that backreaction does not necessarily force us to a small value for \(3q\. 
If M ~ Mpi, then \(3q\ ~ H/M p \ ~ 1CT 6 , while if the new physics is at scales smaller 
than Mpi, we can get a large value for \{3 \. In general we require H < M for inflation 
to be described within the regime of effective field theory, and we must also require 
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that M < Mpi. Nevertheless for the reasonable choice M ~ 10~ 4 M p i we can obtain 
\Po\ ~ 1 for e,r)' ~ 10~ 2 . There are bounds on how large can be coming from direct 
observations of the power spectrum, but for sufficiently small or weakly fc-dependent 
Pk we can easily evade these [29]. 

To summarize: we treat inflaton fluctuations as being described by an effective 
field theory valid at momenta and energies below a scale M. This forces us to the 
notion of an earliest time 770 at which we can trust this theory. We then argue that 
a more natural set of states to use to describe inflaton fluctuations are Bogoliubov 
transforms of the BD state and we ensure that they remain Hadamard by cutting off 
the Bogoliubov coefficient /3& for k > Ma(r) ). These coefficients are further constrained 
by the requirement that the backreaction on the inflating background is negligible. Our 
assumption that the states are Hadamard guarantees renormalizability (in the effective 
field theory sense). We now turn to the calculation of the bi-spectrum using these 
states. 

3. The Inflationary Three Point Function 

Our goal in this section is to show what impact the modifications to the quantum state 
of inflaton fluctuations described in Sec. 2 have on the momentum space (as well as 
£-space) structure of the three point function. To study fluctuations about the FRW 
background, we follow Maldacena [7] and use the ADM foliation of the spacetime as 



where N is the lapse function and N l the shift function. The gauge invariance of 
the Einstein action needs to be fixed and this can be done in one of several ways. For 
instance we can set the fluctuations in the inflaton field to zero 5$ = and parametrize 
the metric fluctuations as: 



The curvature fluctuation ( is useful since it is a gauge invariant quantity and remains 
constant outside the horizon. Note that this definition of ( differs by a sign from 
that often found in the literature and this is relevant to the sign of the three point 
function since the bounds on f NL are not symmetric [30]. For superhorizon wavelengths, 
gradients in ( can be neglected so that the ( can be viewed as shift in the time variable. 
Maldacena's calculation involved using this gauge and parametrization then expanding 
the Einstein action to cubic order in (. 



ds 2 = -N 2 dt 2 + hij (dx i + Wdt) (dx j + N j dt) , 



(3.1) 



= exp (2Ht + 2() h ih det h = 1. 



(3.2) 
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Another gauge that can be quite useful moves all the scalar fluctuations into the 
inflaton field: 

$ (x , t)=<j> it) + 50 (£, t) , hy = exp (2Ht) hij. (3.3) 

In this gauge, we compute the bi-spectrum by first calculating the three point function 
of the fluctuations 5(f), (5(f>8(f>8(f>) and then relating this to (CCC)j it is this latter quantity 
that encodes the observable non-gaussinities in the CMB. To go from 50 to ( we have 
to take into account the non-linear evolution of 50 outside the horizon which involves 
going to second order in the perturbations. This gives rise to a non-linear term of the 
form (5050) (5050) to (CCC) which serves to cancel the time dependence outside the 
horizon coming from the inflaton fluctuations so that ( remains time independent. To 
lowest order in slow roll parameters, which suffices for our needs, the relation between 
50 and ( is 

C = -?. S <l> + 0(e, V )(5-6A , (3.4) 





where e, rj are the slow roll parameters. We will use this latter gauge for our calculations, 
and only compute the leading order behavior in both the slow roll parameters, as well 
as in H/M. This means in particular that the fluctuation modes will be those defined 
on de Sitter. 

As discussed in Sec. 2, in a cosmological setting, the relevant expectation values 
are of the in-in type, as opposed to the in-out ones used to compute S-matrix elements. 
There is a formalism in place for doing this [23] and it has been further elaborated in 
Refs. [11, 24, 31, 32]. We use these techniques here to compute the three point function. 

In the canonical version of this formalism used by Maldacena, we first construct 
the interacting Hamiltonian, Hj, in the usual manner. The equal time tree level con- 
tribution to the three point function is then given by 

(C(zi)C(z2)C(z3)> = - % f <V <[C(*i)C(z2)C(*3), W)]>o; (3.5) 

using the reality condition of the equal time product this is the same as 

(C(*i)C(z2)C(a*)> = "2^e ^V*<C(*i)C(*2)C(* 3 )/W))o) ■ (3.6) 

The lower limit on the integrals corresponds to the beginning of inflation i]q. This is 
not an artificial cutoff, but rather it reflects our choice of a gaussian initial state. In 
practice, this cutoff is crucial if we are to make sense of the integrals when initial states 
other than the Bunch-Davies one are chosen. 
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The free field correlators (((xi)((x2)((x3)Hi(r]')) can be computed via Wick's 
theorem, where contractions are replaced by Wightman functions 4 



(O 1 C(x 1 )O 2 C(x 2 )O 3 ) - (C(zi)C(z 2 ))o(Oi(9203)o + . . . . (3.7) 

Tadpoles can be removed by normal ordering the interacting Hamiltonian i.e. in prac- 
tice neglecting self contractions. The choice of vacuum is then equivalent to the choice 
of Wightman functions 

/r\ 3 k 

j^y 3 VMVMe^-^. (3.8) 

The modes V k satisfy the free field equation 

V'l + 2-V k + k 2 V k = 0, (3.9) 

z 

where z = a<p/H. We choose them to be normalised as 

V k - -i=- (a k e-^ + (3 k e ik *) , (3.10) 
y2kz 



in the limit \ki]\ ^> 1. The three point function will then be given by a causal time 
integral over products of Wightman functions and their derivatives. We are now ready 
to turn to the specific calculations. 

3.1 Scalar Minimally Coupled to Gravity 

Consider now the situation treated by Maldacena, where the inflaton is a scalar field 
minimally coupled to gravity. The result contains two contributions, one from a local 
(in time) nonlinear field redefinition, and the second from an integral of the form in 
Eq. (3.5). Since local redefinitions will not contain the enhancement terms we are 
interested in, we can concentrate on the non-local contribution. Maldacena [7] showed 
that up to these redefinitions the action contains the following three point interaction 
for ( (we have reexpressed his result in conformal time) 



S {3) = Jd V d 3 x a 3 M- 2 U) HC 2 d~X, (3.11) 



where denotes the cosmic time derivative of the inflaton zero mode. We use it 
instead of the conformal time derivative since the slow-roll conditions are easier to 



'This is the "in-in" version of the normal Wick's theorem for time ordered products. 



- 12 - 



express in terms of 0; ~ y2e M p \H. The associated interacting Hamiltonian is 

Hi = — J d 3 x a 3 H^' 2 d~ 2 ^'. Note that this conversion is not entirely trivial 

due to the time-derivative dependence of the interaction. Substituting in the general 
expression Eq. (3.6) we obtain 



4 



(CfcCfeCfe) = -*(27r) 3 5 3 5>) M^H^dr, a 3 ( V ) ^ d v Gl(0, V )dr,Gl (0, V )d v G> 3 (0, v ) 

+ permutations + c.c. . (3-12) 

In writing the above expression, we have evaluated the three point function at rj = 0, 
i.e. when the modes are well outside the Hubble-horizon. Here we see explicitly how the 
product of Wightman functions appears. For the BD vacuum the Wightman function 
is given be 

GI(V,V') = ^^{l + ik V ){l - ik^e-^r (3.13) 
For the first argument taken to be well after horizon crossing, we have 

and so (in the limit r] — > — oo) 

(C S ,C fe C Sa > = fc^-gMi, ( 3,5) 

where k t = k\ + k 2 + k%, and fcj = \ki\. 

Suppose we now modify the mode functions: Uk — > Vk- We get two types of correc- 
tions: one from the modifications to the norm of the positive frequency modes akUk t 
which just changes the overall normalisation and shape dependences by an amount 
similar to that for the power spectrum. The second type of corrections which are sup- 
pressed by powers of j3k are more interesting. One such correction will again modify 
the over amplitude, but the more relevant corrections are those that change the argu- 
ments of the oscillatory exponentials. For instance, to linear order in (3% we have the 
correction 

where we have defined kj — k t — 2k j. Performing the integral we have 

E_ o ff6 k 2 k 2 k 2 / - \ 

ki) - . Y ( 1 - e ik ^) + c.c. . 

(3.17) 
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What makes this contribution interesting are the denominators proportional to kj 
which appear. In the standard calculation in Ref. [7], the denominator only involved 
k t = ki + k 2 + k 3 which could never vanish unless all the /c's vanished, which would 
then force the numerator to vanish as well. However, the mixing between positive and 
negative energy BD states in our modes now allows for vanishing denominators (kj = 0) 
for certain non-zero values of the momenta, which enhances these contributions relative 
to those for the BD vacuum. In practice there is no divergence since the exponential 
factor is unity in the kj — > limit. However, this cutoff is only relevant for kj ~ l/i]o- 
Thus if we evaluate Eq. (3.17) for these special triangles where kj = 0, we get 

AfefeCj,)^. = (arWE k) n( ^ ?) Jfj, ^ MA>, (3-18) 
and so relative to the standard result we get an enhancement factor 



(CjfcjCfcoCfc,) 



Vgp *> i*" fc "«i = ^^m- < 3 - 19 > 



We have implicitly assumed that Im/^ ~ \flk\- Note then that although the result is 
suppressed by a factor of /3k, this is multiplied by the ratio of the physical momentum at 
the beginning of inflation to the Hubble scale, which may be as large as M/H (at least 
for those modes which have been excited). In practice however, we cannot measure 
the full three point function but only the two dimensional projection of it encoded into 
the CMB. As a result we must effectively average this result over angles which we will 
do in Sec. 4. On doing so we essentially lose a factor of \krj \ and so this result is 
not enhanced but will be of the same order as similar small corrections from the local 
redefinitions and to the two-point function. 

Nevertheless, for interactions which scale with larger powers oil/ a the enhance- 
ment effect can be significant and we shall discuss these in the next section. The 
enhancement will occur when e.g. k\ + k 2 = k%. From the momentum delta function, 
we also have that \k\ + k 2 \ = k 3 so that only "flattened" triangles with two of the 
vectors being collinear will be enhanced. Whether or not the effect is observable then 
has to do with whether these triangles have sufficient "measure" when the bi-spectrum 
is converted to £-space. We will discuss this further below. 

3.2 Higher derivative interactions 

Now let us turn to the model of higher derivative interactions discussed in [6]. It 
contains the following dimension 8 correction to the effective action for the inflaton: 

A 



Cl = ^sW^ )2)2 - (3 - 20) 
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This type of correction is regularly considered in modifications to the kinetic term of 
the inflaton such as in k-inflation [33], or DBI inflation [34]. Since we do not want to 
be wedded to a specific scenario, for the purposes of the present discussion we shall 
concentrate on standard slow-roll inflation; however, it is straightforward to extend the 
arguments to these more general cases. 

Before looking at its contribution to the three point function, we shall make some 
cautionary remarks about the tree level contribution to the two point function coming 
from this interaction. Expanding $ = (f)(i]) + 5(f) to second order in 5(f) and then 
converting from 5(f> to (, this term contributes the following term to the Hamiltonian 

H?^ = ^ J (3^,2 _ m)2 j 

J dW-^(3(C' + arC) 2 -(^C) 2 ), (3.21) 



A 



4M 4 

where r = dln ^ H ^ ■ Since r is both slow roll suppressed and its contribution to the 
action subdominant at subhorizon scales we can ignore it for the present discussion. 
Since this term is quadratic in the fields, we could, in principle, simply incorporate it 
into the free Hamiltonian. On the other hand, we could also treat this as an interaction, 
just as for the three point function. Doing this would yield the following correction to 
the C equal time two point function evaluated at late times 

A(C(*i)C(* 2 )> = -^fj^Wjp (3(d v G>(0,v)) 2 -k 2 G>(0,v) 2 )+cc. (3.22) 

On considering the terms linear in (3 k , this appears to give similar enhancements as 
described in the previous section. In this case the enhancements will not get washed 
out by averaging over angles. However, this is clearly a fake since had we incorporated 
this term into the free Hamiltonian we would not have seen any such affect. The 
resolution is that this term can be removed by a renormalization of the initial vacuum 
choice flk — > Pk + X5j3k- 

To see where the problem arises we can always resum two point vertices by ab- 
sorbing them into the definition of the free Hamiltonian. On doing so, the Wightman 
function will be constructed from free modes satisfying an equation of the form 

V£ + 2-V k + c 2 s k 2 V k = 0, (3.23) 

z 



where 



2 (1 - A0 2 /(2M 2 )) 



ci = ± ^ '- L . (3.24) 

s (1-A0 2 /(3M 2 )) V ; 
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Here we see the modified sound speed characteristic of models such as k-inflation [33] 
and DBI inflation [34] that make use of these higher derivative operators. At subhorizon 
scales the WKB approximation is valid, and so can write the general solution as 



= /7 r— r (a k exp(-ic s \k\ri) + (3 k exp(+ic s \k\r])) . (3.25) 
y/2c s kz 

We can now see the reason for the apparently large corrections to the power spectrum. 
If we expand these modes out to 0(A) we get 

exp(-tc s \k\ V ) = exp(-t\k\ V )(l - i\k\v^ + 0(A 2 )), (3.26) 

which is a poor expansion for sufficiently large \kt]\. From now on we shall work with 
the resummed mode functions V k which are well behaved at subhorizon scales, which 
has the effect of removing the correction in Eq. (3.22); since the correction to the three 
point function computed below is already O(X), we can omit these corrections to the 
modes. 

Having taking care of these issues, we can now turn to the calculation of the three 
point function in the presence of these higher derivative interactions. Expanding the 
interaction Lagrangian to third order in the fluctuations S(f> and then converting to (, 
we find 

Hf = - J 6?xa^^ - OM?) , 

where again we have neglected the r term. Note that although the original operator 
was dimension 8, this operator is really dimension 6 in terms of 5(f>, but it is additionally 
suppressed by y/e (HM p i/M 4 ). 

We can now compute the corrections of interest to the three point function. We 
compute to lowest order in A, H/M, and e to find 

(CfciCteCte) = (27r) 3 5( 3 ) A (*1>*2,* 3 ) 

A (\ k 2 , k 3 ) = J" d V a {d v Gl(0, V )d v G> 2 (0, v)d v G> 3 (0, V ) + 

+ (h ■ k 2 )(G>(0,rj)G> 2 (0, ^^(O, rj)) + perms.) +c.c. (3.28) 



- 16 - 



The corrections we want are again those linear in (3k- These give rise to the following 
contribution to the three point function ( k . 



E k j n (2fc 3) 02M4 E ^ l o dve^Sjih, k 2 , k 3 , r,)+cc, 

(3.29) 



where 



Sj(ki, k 2 , fc 3 ,77) = -k t [Uiki) + (3.30) 
(-77) (tf(k j+1 + k j+2 ) + k j+1 k j+2 (k j+1 - A; j+2 ) 2 (A; i+ i + £; j+2 ) - tf(k 2 j+l + fc? +2 ) 
-k 2 j(k j+1 + k j+2 ){k 2 j+1 + k j+l k j+2 + fcf +2 )) + (-77)% (n^) (A; 2 - 4k j+1 k j+2 ) , 

where j is defined modulo 3. As in Sec. 3.1, we see that these results are enhanced 
for the flattened triangles where two of the vectors are collinear so that fc, vanishes for 
some i. Note that the first and third terms in the 77 expansion of Sj vanish and so the 
dominant contribution evaluated on a give flattened triangle kj = is 

A<c £ ,< £ ,c S3 >|,. = - W*" (es) ukn^^ < 3 - 31 > 

X 7T {4kj+ik j+2 (k j+1 + k j+2 )(k 2 +1 + k 2 +2 + ^+1^4-2)) • 
For these specific triangles the enhancement factor relative to the BD contribution is 

A(C 3 



(C 3 



~ l&INol, (3.32) 

flattened 



where we have assumed TZe((3k) ~ \flk\- Again as will be explained in the next section 
on going to I space we effectively lose one factor of \krj \, but one enhancement factor 
remains. Thus we can give an order of magnitude estimate for f^L for these triangles 
to be 

fNL I flattened ~ T7±\Pk\ [ „(„ \ tt • ( 3 ' 33 ) 



M 4 \a(7 l0 )H / 

For backreaction to be under control, the largest reasonable value for \/3k \ is y/erfHM p \/M 2 
while for the effective field theory to be valid the largest value allowed for k/a(r) ) is 
M. Thus the maximum expected contribution to f^L is 

/^flattened ~ vV^^ [h ) ^ ^ [m^J {~M ) ' ^ 
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We see then that for reasonable values of these parameters e, 77' ~ 10 2 , H/M v \ 
1CT 6 we get 



6.3xlO- 4 ^J , (3.35) 



which is in the range 1 to 100 for M in the range 6.3 x 10~ 4 M pl to 2.5 x 10" 4 M pl . In 
pushing into the limit /nl ~ 100 we have \(3k\ ~ 1-6 x 10 _1 which may or may not be 
already observable in the power spectrum depending on the precise k dependence of 
the Pk- In practice these bounds are overly restrictive since we know that models such 
as DBI inflation do not have to satisfy the usual slow roll restrictions, and in fact the 
more general case was considered in Ref. [17]. In these models the usual contribution 
to the non-gaussianity may be large, and so this specific contribution will be further 
enhanced. In general the enhancement factor for these specific triangles relative to the 
usual contribution is bounded by 

1 100y/\fa\. (3.36) 



I'nl 


flattened 


Jnl\ 


usual 



From an observational point of view, one might think that since only specific triangles 
give the enhancement, this effect might be buried in the noise. However, consider the 
I = 500 total modes which have been observed by WMAP with signal to noise greater 
than 1. In total, the bi-spectrum is made of (500) 3 points (modulo symmetry factors). 
Of these, the total number that satisfy the triangle inequality are roughly (500) 2 and so 
are down by a factor of 500. This would then imply that the signal to noise coming from 
these will thus be down by a factor of 1/ V500. Thus, roughly 5 percent of the signal 
will come from these modes, and so to a crude approximation ~ 10 2 a/|/5a;| enhancement 
will effectively only amount to a ~ IOa/I/^I enhancement. In practice it is necessary to 
reanalyse the data with an appropriate template along the lines of Ref. [30] , and this 
will improve the constraints on this contribution. 

3.3 Backreaction from higher derivative interactions 

We have already given a simple estimate for the absence of backreaction based on 
computing the expectation value of the free field stress energy. However, as soon as we 
choose a non-standard initial state, and a non-zero gaussianity develops we will also 
generate additional cubic and higher order contributions to the stress energy. As usual 
we can infer these via the in-in formalism [23] so that the expectation value of the stress 
energy at any time 77 is given be 

(T^fax)) = (Te i % d * H 'T lu ,(T,,x)Te- i %> drfHl ) , (3.37) 
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which gives to first order in the interaction 

(7^(77, x)) = -2TZe (i J V dj/^fa, x)H I ) ^j , (3.38) 

where it is important to not that the T^ v inserted on the right hand side is defined 
in the interacting representation. Thus in particular a contribution to T^ v cubic in 
the interaction muliplied by the cubic vertex in the Hamiltonian will give a non-zero 
expectation value. If we choose the Bunch-Davies vacuum then for the usual reasons we 
can safely ignore this backreaction contribution since the vacuum is de Sitter invariant 
and it will at most contribute to a renormalization of the zero point of the potential. 
If, however, we choose instead an excited yet still gaussian initial state, the early non- 
gaussianity that is generated may cause a problem. 

Adding to the action the interaction g^j((90) 4 gives the following contribution to 
the stress energy (in the Heisenberg rep.) 

AT£ = -^WfW + ^W) V- ( 3 - 39 ) 

In particular, expanding to cubic order around an inflating solution 4>{rj) (and neglecting 
backreaction of the metric which is a reasonable approximation at sub-Hubble scales) 
we get the following cubic contribution to the energy density 

A (3 )P" = 2]^s^ ( 3 ^' 2 - W) 2 ) • (3-40) 

On converting to the interacting representation we unsurprisingly obtain the interacting 
Hamiltonian density 5 

A( 3 )P = jj^s^' (V 2 - W) 2 ) • (3-41) 

In any case rather than performing the full calculation we may simply note that at 
sub-Hubble scales 8(j)' and di5<f) scale the same way so the effective contributions to the 
pressure and stress energy take the form (where we have defined x = a<50) 6 

\2l2 rv r 1 

A <P I j../ / j3„/ 1 /. /3/ \. /3, 



' '/() 



A <» p 53 imf) L d "7 ^'^"(vM'.^ (3-42) 



M 8 a 6 (r]) 



J'drj'J d'x'-^id^Hix'^-rj')) 3 , (3.43) 



5 To get from one to the other one first writes the Heisenberg expression in terms of <j> and it, and 
then replace it by its free field value it = a 3 (j>. 

6 In fact the final answer is smaller since the Lorentz invariant combination <50' 2 — (diScf)) 2 gives an 
additional H/M suppression. 
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where 

d 3 A; 1 



+ al(3 k S k ^+V) + a k (3le- m{r > +rf) + |/3 fe | V 1 * 1 ^^) . (3.44) 
This gives at worst a contribution of order 

A, n ~ A2( ^ 2 7?r r dn 1 —^— f f frjjk _i_ jkJfl* p -2ifeir? -i(-fci+fc 2 +|fei+fe 2 |)(W) 

(3.45) 

to first order in As explained in Sec. 2 the presence of the oscillatory factor e ±2lkiri 
washes out the early time contribution to this integral. The leading contribution is 
then the O(Pl) contribution which is bounded by 

Ap (3 ) < |A>| 2 2 - (3.46) 

To obtain this estimate we can simply count the positive powers of k and rj in the inte- 
gral and replace them with their maximum possible values a(r)o)M and rjo respectively, 
and then dividing by one factor of \krj \ from the angular integration. Note that since 
we also require 2 < M 4 for the validity of the effective field theory this contribution 
is necessarily smaller than the free field contribution. 

Similarly at quartic order we have in the Heisenberg representation 

A (4)P H = ( 3( ^ _ 2 (50') 2 (W) 2 _ . (3.47) 

This gives a contribution of the form 

A ( 4,P - M^£*ff dV^(x' 4 (^) X 'V,*')>o (3.48) 
jT dr/ J dV-^y (d^Hix', 77 - i)f . (3.49) 

As before the dominant contribution is bounded by 

AP(4) < \Po\ 2 M\ (3.50) 

According to the general argument given in Sec. 2 higher order contributions will never 
give a contribution larger than this as long as we remain in the validity of the effective 
theory. 



A 2 

M 8 a 8 (7]) 
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4. Measure issues and the 1-space bi-spectrum 



The results of the previous section show that for flattened triangles the three point 
function can be considerably enhanced relative to the usual case. However, we do 
not actually measure the three dimensional three point function directly in the CMB. 
Rather, we measure its two dimensional projection as encoded in the a; m 's. 

The relevant question then is to what extent these enhancements can be seen in 
the CMB 1-space three point function. Once the temperature fluctuations have been 
decomposed into spherical harmonics, 

AT 

-f-fr) = ^2 a imYi m {n), (4.1) 

lm 

we can construct the angular averaged bi-spectrum [35] 

B(h,l2,h)= ( 1 2 3 ) {o-hrn 1 0-l 2 rn 2 a hrnz) ■ (4.2) 

\ mi ^2 rri3 J 

mi,m2,m3 x ' 

Expressing the three point function for ( k as 

(C&C&Cfe) = (27r) s *(£ h) A(h, k 2 , k 3 ), (4.3) 

i 

the three point function of the a\ m is given by (in this section we will make use of the 
unit vectors defined via ki — |/uj|fij). 

/A3U r] 3 kn H 3 ?q 
(27T)3 (27T)3 (27T)3 Krn^)Krn^) Y llrn 3 (n 3 ) 

(27r) 3 5( 3 )(^^)^(fc 1 ,fc 2 ,fc 3 )A^(A; 1 )A^A; 2 )Aj(A;3), (4.4) 

i 

where Af(k) are the radiation transfer functions. We are interested in a specific com- 
bination to A that comes from flattened triangles. The contributions to the three point 
function linear in f3k can be written in the form 

A(h,k 2 ,h) = $>4(^r +1 f°d n (-riye^V(k j ,k j+1 ,k j+2 ) + c.c., (4.5) 

3 Jr >° 

where j is defined modulo 3. Let us consider the properties of the integral 

r° 

I(k t ) = I dr](~r]) s e ik ^. (4.6) 
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For \kiT]o | ^> 1 the oscillatory nature of the integrand heavily damps the power of f] 
and as a result the integral scales as (l/iki) s+1 . On the other hand for l&^ol <^ 1 
the exponential is irrelevant and the integral gives I (hi) = (— r] ) s+1 /(s + 1). We are 
interested in effects coming from sub-Horizon scales where \kiT] \ ^> 1, and so the 
dominant contribution to the integral will come from the regime where the triangles 
are flattened and ki ~ so that the integral is given by its maximum value I{k,j) = 
(— Vo) s+1 / (s+ 1)- However, this will only occur for a special range of angles, essentially 
for 

n t .n 3 + 1 « ( \~ fcj ° 2 . (4.7) 

To perform the integrals, let us first integrate out kj + 2 via the delta function, and 
then integrate over hi using n i+ i as a reference axis. On doing so we get 



(a hmi a hm2 a hm3 ) = ^?i h+h+h J ^ / / J ^ 

^(n*)^ +1 m J+1 (^ \kj - k j+1 \) x 

A?(^)AJ +1 (% +1 )AJ +2 (|^ - k j+1 \) + . . . , (4.8) 

where the ellipsis represents the terms arising from the complex conjugate in Eq. (4.5). 
At this point we can make the approximation that since the integral I(kj) is so strongly 
peaked at kj = we may replace the product of spherical harmonics with their values 
for the flattened triangles set by kj — 0, namely 

^■mj (^i+l)^* +2 m j+2 (^i+ 2 ) ~^ ^*m 3 ( — ™i+l)^* + im J+1 (%+l)^* +2 m j+2 

(4.9) 

Similarly we may replace 

AJ +2 (|% - k j+1 \) - AJ +2 (|A;, - fc i+1 |). (4.10) 
In this approximation, the integral I{kj) gets directly integrated over angles 



/ 



&%I(k 3 ) = 2tt f d V (- V y—l— 2 (e 2ik ^ (1 - i(kj + k j+1 ) V ) - (1 - i(kj - k j+1 )rj)) , 



(4.11) 

where we have assumed kj > kj+i as is necessary for the existence of this triangle. In 
the limit \k s r) Q \ ^> 1 this integral is given approximately by 

d 2 V(M - M t'7* J " +l) HE.)'- (4-12) 
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So we have 

-2 



^; mj (-n J+ i)^ +im3+l (n j+1 )^; +2m , +2 (n i+1 )^ 2 ^+i - %) H w x 

SKjK j+1 

V(k h k J+1 , \kj - kj^Alik^AlJkj^AlJlk, - k j+1 \) + .... (4.13) 
Using the fact that Yi m (—n) = (—l) l Yi m (n) we get 



{ai imi a hm2 ai 3m3 ) ~ 2_^l 47r J * I -U y ^ 27r )3 y fc fe (2^)3 P k i y h,. 
2irk t (k j+1 - fcj) 



i-iktmyvfa, k j+1 , \k 3 - k^DAi^AKkj+JAf^kj - k j+1 \) + . . ., 



where G™]'^'™ 3 is the Gaunt integral [35]. Finally then, these contributions to the 
bi-spectrum take the form 

B{hMM) = (4.14) 
^ / (2/ 1 + l)(2/ 2 + l) (2/ 3 + l) /<i / 2 h\ ( ,, 3 , 1+h+l3( _^ h f dkjk] f dk J+1 k] +1 
\ V 4vr lo O) W { 1} J (2vr)3 J kj+i<kj (2*)s 

^ /^ ( ^ + ~ fej) (-iMo) s P(%, I A* - % + i|)AJ(%)AJ +i (A; j+1 )A^ 2 (|A;, - fc i+1 |) + c.c. 

The crucial point to note is that whilst the original integral naively scales as \ktT]o\ s+1 , 
the above bi-spectrum only contains a factor \kti]o\ s . In short, in performing the 2-d 
projection we essentially have to smooth the triangles over some finite resolution of solid 
angles. This softens the size of the enhancement by one power of \kr] \ = k/(a(r)o)H). 
Nevertheless as we have described in Sec. 3 the remaining effect can feasibly be larger 
that the usual one. The radiation transfer functions are dominated by the multipole 
I ~ kd,Lss where cIlss is the distance to the last scattering surface, and so it is clear 
that the product Af (kj) Af (kj+i) Af. +2 (\kj — kj + ±\) is dominated by the multipoles 
which saturate the triangle inequality, namely lj = + lj+2- The precise form of the 
bi-spectrum will depend on our model for the k dependence of the (3^, of which we can 
only guess the form, but the characteristic peak for the flattened triangles should be 
sufficient to distinguish this effect from the usual contributions to non-gaussianity. Note 
that in the squeezed limit this leading order contribution vanishes, and a subleading 
term takes over which is less enhanced. Thus this contribution to the three point 
function should be clearly distinguishable from the squeezed triangles contribution. 
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5. Higher irrelevant operators, N-point functions 



Now let us consider higher order operators and their contributions to various N-point 
functions. It is possible to make fairly general statements based on the scaling with 1 ja 
as to the magnitude of the terms. For the three point function let us first consider the 
effects of higher order operators of the form M 2™+4 (( V0) 2 ) n . Naively since these are 
increasingly irrelevant operators we might expect them to scale with increasing powers 
of 1/a. However, expanding to cubic order and defining x — a ^(p (this is usual since x 
is pure oscillatory at sub-Hubble scales with no scale factor dependence) we find 

2n-3 



3-pt u 1V1 




^ (cix' 3 + c 2X '(V X ) 2 ). (5.1) 



So all the higher order operators will have the same scaling with 1/a and are just 
additionally suppressed by <fi/M 2 . In particular for DBI-inflation and k-inflation we 
will still have an interaction scaling as 1/a 2 and so similar enhancements to those 
described. Similar statements can be made for the N point functions where we find 

/ . \ 2n-N 

(W) 2 )'V P . - U-A (C3X'"+C4X" V - 2 (VX) 2 + ...). 



(5.2) 

To increase powers of 1/a we must go to genuinely higher derivative interactions (rather 
than powers of first derivatives) such as (V a Vb0V a V b 0) 2 . Note that terms con- 
taining D(p can be neglected as they are redundant couplings which can be removed 
by local field redefinitions. Although the correlation functions are not invariant under 
these redefinitions, it is precisely because they are local that they will not give rise to 
any interesting effects at late times. Calculating the three point function we have 



(V a V 6 0V«V fc 



! 4>(d a d bX ) 2 x" « T^-Mdad b x)%d*x- (5.3) 



3 _ pt M 8 a 5 " v " * M 8 a 



So we see that this term scales as 1/a 5 . Similarly, there are yet higher derivative inter- 
actions that we can construct that scale as higher powers of 1/a. However, these naive 
scalings may not directly represent the enhancements that arise due to the cancellation 
of numerators evaluated on the flattened triangles. Further more these higher order 
derivative terms are necessarily additionally H/M suppressed which typically cancels 
out the extra 1/a enhancement . Nevertheless the general form of the enhancement ef- 
fects for the three point function implied by effective field theory will take the schematic 
form (to first order in (3k) 



v™{^%>t) +o W' (5 - 4) 
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where the function / is a dimensionless function which has a Taylor expansion in its 
first argument and similar relations holding for the N-pt functions. Thus we see that 
if f3 k 7^ 0, the validity of effective field theory requires that k/(Ma(r) )) <C 1 (for the 
Taylor expansion to be valid) as we stressed from the outset. With sufficient e-folds 
this can be violated and so we discover directly the transplankian problem [26], namely 
the statement that effective field theory does breaks down at the beginning of inflation 
if the initial vacuum state is chosen to be nontrivial. However if j3k = there is no 
reason to doubt effective field theory even when k > Ma. The physical reason is clear; 
the BD vacuum contains no particles, so that there is nothing that physically carries 
transplankian energies, and so no transplankian effects should be expected. 



6. Conclusions 

Our main conclusions can be summarized as follows: Choosing an excited gaussian 
(and Hadamard) initial state at the beginning of inflation can give rise to additional 
contributions to the bi-spectrum from sub-Hubble interactions. They are significant 
for flattened triangles k\ = k 2 + k 3 (and permutations), which in / space saturate the 
triangle inequality l± = l 2 + h (and permutations). Schematically the corrections from 
a given operator calculated in an excited state, relative to the equivalent BD state 
contribution take the form 



A(C 3 ) 



(C 3 ) 



~\MkVo\ s+ \ (6.1) 

flattened 



where the integer s depends on the precise scaling of the interaction with inverse scale 
factor. The specific contribution to /jvl for higher derivative interactions is of the order 



A(C 3 ) 



(C 3 



flattened 



The requirement that the excited inflaton quanta do not spoil the slow roll evolution 
imposes the bound \(3k\ < \ferf while the validity of effective field theory bounds 
< M. Thus the largest expected enhancement factor is 



A(C 3 ) 



(C 3 



j—HMj* My ^M pl 

flattened " ^ ^ \h) = ^ If ^ 



which can easily be of order 10 4 . Whilst it may be possible to observe this enhancement 
in a genuinely three dimensional probe, the CMB bi-spectrum is only sensitive to a two 
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dimensional projection of this. On converting to I space we must average over angles 
which gives a reduction of one factor of \kr] \ so that 



A/. 



NL 



~ \MkVo\ s - (6.4) 

flattened 



fNL 

In particular for the higher derivative interactions, the largest expected value is 

A/^Uene d ~ ^^f^ ~ ~ WW* ■ (6.5) 

In this particular case we need \f3k\ > 2.5 x 1CT 2 to obtain /jvl > 1- Whether this 
departure from BD can already be detected in the power spectrum, depends on the 
precise k dependence of \/3k\- Note that in the extreme case in which (3 k is constant over 
the range of scales probed by the CMB, the only affect of f3 k in the power spectrum 
would be to renormalize its amplitude. On the other hand a stronger /c-dependence 
could be probed in any spectral tilt present. Nevertheless to get an enhancement 
relative to the usual result we only need \/3k\ > 10~ 4 and so clearly models which 
already give a large non-gaussian contribution will be the most interesting in terms of 
seeing the enhancement. 

It is clear from our discussion that typically these effects will only be seen if the 
number of e-folds of inflation is low, i.e. little more than the require ~ 55 — 60, since 
we require that the modes corresponding to the observable scales today had physical 
momenta less that M < M p \ at the onset of inflation. Thus the effect is in the exact 
opposite region to the transplanckian problem which typically requires ~ 70 e-folds 
for the physical momenta to reach M p \. We can imagine a two-stage inflation scenario 
however where inflation occurs before this period for a possibly large number of e-folds 
with the usual Bunch-Davies vacuum, and then some non-adiabatic process excites the 
vacuum state up to scales set by the Hubble constant at that time ~ M (e.g. a sudden 
change in the trajectory of the inflaton in a multidimensional landscape). Then if a 
second period of lower scale inflation occurred during which the fluctuations relevant 
for observable scales today were created, these would arise from this excited state and 
show the signatures we have described. 

The essential feature of our result was anticipated independently in Ref. [17] where 
a general calculation was performed for theories of the form p(X,<f>). However these 
authors did not give any estimate as to the possible magnitude of the effect which 
we have here emphasized is strongly constrained by backreaction. Furthermore as 
explained there is little chance of observing the full enhancement for the flattened 
triangles, but rather only smoothed one which is further suppressed by H/M. When 
dealing with theories of the form p(X, 0) there is an implicit assumption that one can 
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take seriously the effective field theory beyond the normal expected range of validity, 
by taking seriously the effect of higher derivatives. By contrast in this article our 
perspective has been that even in standard slow roll minimally coupled inflation, these 
corrections will exist due to standard effective field theory arguments, and we have 
remained strictly in a regime in which they are small so that standard effective field 
theory is valid. To properly analyze the case considered in Ref. [17] it will be necessary 
to redo the backreaction analysis [37]. 

Irrelevant operators may have very relevant physical consequences, and indeed as 
we see here, because the relevant contribution to the three point function is naturally 
suppressed, the first irrelevant higher derivative operator can give a larger correction. 
For these effects to be observable the cutoff scale of the theory must be less than 
the Planck scale. Nevertheless, an effect can easily be present for reasonable scales 
H ~ 10 13 ~ 14 GeV and M ~ 10 15 ~ 16 GeV. We may note that super-horizon approaches 
to the calculation of non-gaussianities [36] will not be able to describe the effects we 
have considered precisely because they arise from sub-Hubble scale interactions. 

From a theoretical point of view, we can use this calculation to get a better un- 
derstanding of the special nature of the BD vacuum as the true vacuum of interacting 
quantum fields in de Sitter. What we have shown here can be reinterpreted as the 
statement that the BD vacuum is self-consistent in that if we choose it to define a 
Gaussian initial state, although a non-zero three-point function will be generated, it 
will be highly suppressed both by slow- roll parameters as well as by factors of H/M. 
Any initial deviations from the BD state will induce potentially large corrections to 
the three-point function (although it is important to note that these are still under 
perturbative control). Keeping the three point function small enough to be consistent 
with observations then requires a fine-tuning of the initial state (via the /V s )- 

We believe these results also offer a new window to address the transplankian 
problem, since for the reasons we have explained, higher N-pt functions can be a better 
probe of the transplankian era that the two point function, provided at least that we 
start out in an initially excited state. On the other hand it is likely that if the initial 
state is taken to be the equivalent of the BD vacuum in the transplankian region, these 
effects will be suppressed at least by H/M. 

In conclusion then we see that the CMB bi-spectrum will be a most useful tool for 
probing the inflaton initial state, and can contain richer information than the power 
spectrum alone. Initial state effects will be present in flattened triangles, and a direct 
observation for an enhancement in these triangles would provide the strongest window 
into the physics at the beginning of inflation. 
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